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Abstract 


This  technical  treatise  has  exemplified  a  response  spectrum  method  that  is  particularly  suited 
for  spring-mass-damping  systems  excited  by  one  or  more  base  or  ground  disturbances 
characterized  by  displacements.  It  offers  an  alternative  to  the  conventional  response  spectrum 
method,  which  is  best  suited  for  a  single  prescribed  ground  acceleration.  This  methodology  has 
also  introduced  a  method  of  how  to  incorporate  hysteretic  damping  through  complex  spring 
stiffnesses.  When  the  formulation  is  expressed  via  base  displacements,  a  simpler  application  of 
the  response  spectrum  method  to  varying  spring-mass  systems  with  attendant  input  excitations 
can  be  realized  as  expressed  as  the  ratio  of  energy  considerations  inherent  in  the  natural  physics 
of  the  phenomena  itself. 
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1.  Introduction 


Hysteretic  damping  is  responsible  for  limiting  the  oscillations  in  mechamcal  systems  such  as  in 
drive  trains  of  modem  tanks.  If  we  look  at  a  linear  one-degree-of-freedom  spring,  mass,  and 
dashpot  arrangement,  the  selection  of  the  proper  values  of  mass,  stif&iess,  and  dandling  constants 
to  produce  an  overdanped  or  critically  damped  system  can  be  effectuated  without  too  much 
difficulty.  If  we  consider  an  n-degree-of-freedom  spring  and  mass  distributed  system  model 
where  Xi(t)  would  be  the  absolute  displacement  and  Y(t)  would  be  the  absolute  ground 
displacement,  then  it  is  possible  for  Y(t)  alone  to  be  specified  such  that  Xi{t)  may  be  determined,  as 
in  the  case  of  the  design  of  a  gas  turbine  engine  so  as  to  resist  shock  and  vibration.  The  usual 
classical  approach  is  simply  to  start  out  with  the  laws  of  motion  with  respect  to  the  base  or 
foundation  firamework  (Beskos  and  Boley  1980).  The  base  then  becomes  fixed  where  each 
respective  identified  mass  will  be  subjected  to  an  inertial  force  niSit)  snd  the  displacement 
variables  become  the  relative  displacements,  viz.,  Z,(f). 

Hence,  we  can  now  proceed  with  a  modal  analysis  solution  using  either  the  acceleration  time 
history,  Y(t),  or  the  acceleration  response  spectram  A(to)  that  is  associated  with  Y(t).  Since  shock 
and  base  line  or  even  ground  motion  do  not  have  repeatable  time  histories  but  do  have  repeatable 
shock  spectra,  then  the  acceleration  response  spectram  method  will  be  of  greater  interest  (Biggs 
1964).  It  is  also  possible  to  formulate  our  original  speculative  problem  with  a  fixed  base  and  with 
only  the  mass  next  to  the  base  subjected  to  a  force  equal  to  knY(f)  and  with  the  displacement 
variables  being  the  absolute  displacements,  Xi(t).  What  is  needed  at  this  point  of  analysis  is  to 
sin5)ly  cast  the  equation  of  motion  with  respect  to  an  observer  at  rest  and  move  the  term  k„Y(t)  to 
the  right-hand  side  of  the  n*  equation. 

The  purpose  of  this  treatise  is  to  follow  through  the  consequences  of  this  latter  formulational 
ideology  such  that  it  can  be  easily  invoked  via  a  suitable  computer  simulation  and  then  modify  the 
formulation  to  allow  the  use  of  the  conqjlex  stififiiess  model  of  hysteretic  damping  for  oscillations 
of  multi-degree-of-fireedom  systems,  especially  inherent  in  the  accessory  gear  box,  reduction  gear 
box,  and  output  shaft  of  a  typical  gas  turbine  engine.  This  type  of  analysis  can  prove  to  be 
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invaluable  for  the  diagnostics  and  prognostics  of  engines  as  illustrated  by  Helfinan,  Dumer,  and 
Hanratty  (1995).  This  work  differs  from  other  investigations  by  invoking  con^lex  stiffiiesses  for 
hysteretic  dancing,  whereas  in  most  analyses,  the  stiffiiesses  are  represented  as  ordinary  spring 
constants. 

2.  Concept  Formulation 


The  systems  considered  herein  are  those  that  can  be  modeled  by  the  matrix  differential 
equation 


MX(t}+KXit}=KJ{t\  (1) 

where  X(t)  is  an  n-dimensional  vector  of  displacements  and  M  and  K  are  nxn  symmetric  matrices 
containing  the  physical  parameters  of  mass  and  stiffiiess  constants.  Let  us  assume  that  the 
eigenvalue  problem  associated  with  (1)  has  a  spectral  matrix  and  a  modal  matrix  such  that 
is  the  i*  eigenvector  with  elements  where  r  =  1,2, . . .,  n.  Furthermore,  also  assume  that 
the  are  normalized  so  that  MAX  ,  =1  for  each  respective  i.  In  order  to  have  orthogonal 

r  ' 

conditions,  we  are  required  to  have 


=  Mi  (2) 

and 

=  Ki  (3) 

where  M,  are  the  modal  masses  and  Kt  are  the  modal  stiffiiesses.  Furthermore,  it  can  be  easily 

shown  that  Clj  =  Ki/Mu  If  we  look  at  the  transformation  of  coordinates  (Brent  1973),  viz., 

X(t)  =  <^q(t),  (4) 
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then  (1)  becomes, 


Mi^q  +  mq  =  KJ. 


(5) 


Then,  if  we  premultiply  by  and  use  (2)  and  (3),  we  then  have, 

M^q  +  K,q  =  ^^K„Y.  (6) 

Also,  it  can  be  shown  that 

=«>„,,,  (7) 

i.e.,  a  vector  con^osed  of  the  n**  element  of  each  Oi.  Then,  (5)  becomes 

+  (8) 

The  response  spectrum  is  the  maximum  response  of  an  oscillator  subjected  to  a  forcing  function 

drawn  from  the  process  of  interest  as  a  function  of  the  frequency  of  the  oscillator.  If  the  response 

of  interest  is  the  absolute  acceleration,  then  the  acceleration  response  spectrum  is  denoted  by 
A(a)),  where  o)  is  the  angular  frequency  of  the  oscillator.  If 

A(&),.)=  (9) 

where  Qi(t)  is  the  solution  to 

M,Q  +  K,Q,  =  K,?,  (10) 

then,  linearity  requires  that 

=  ,  (11) 
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where  qi(t)  is  the  solution  to  (8).  Furthermore,  (4)  can  be  recast  as  (Byme  and  Hall  1973) 


(12) 


i=l 


The  acceleration  of  the  r*  degree-of-freedom  in  the  i*  mode  can  therefore  be  expressed  as 


(13) 


and  using  (1 1),  we  then  have 


(14) 


Also,  if  we  let  (Inman  and  Andry  1980) 


MAX  X 


T,  i 


=  R, 


(15) 


then,  in  going  from  (13)  to  (14),  the  information  about  the  temporal  relationships  among  the 
terms  on  the  left-hand  side  of  (15)  will  be  lost.  This  loss  of  information  requires  that  the  exact 
modal  recombination  procedure,  vk.  (4),  be  suitably  replaced  by  an  approximate  modal 
recombination  rule.  If  we  use  the  square  root  of  the  sum  of  the  squares,  i.e.. 


then  Ar  is  expected  to  be  an  approximate  measure  of  the  maximum  absolute  acceleration  of  mass 
point  Mr  when  all  modes  are  excited  by  input  ¥(()•  The  spring  forces  in  mode  i  can  be  obtained 
from 
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Fr,i=Kr 


(17) 


and,  if  required,  they  can  also  be  combined  by  the  square  root  of  the  sum  of  the  squares  rule. 
Equations  (14)  and  (17)  can  also  be  the  basis  for  a  conqjuter  solution  of  (1)  using  the  acceleration 
response  spectrum  associated  Y(t).  If  A(w)  is  simply  not  specified,  a  separate  program  will  be 
necessary  to  compute  A(ci))  firom  Y{t)  in  accordance  with  (10).  In  the  case  of  multiple,  distinct 
inputs,  as  is  the  case  in  the  structure  of  the  gas  turbine  engine,  it  is  better  to  generalize  (1),  as 
done  similarly  by  Biggs  (1964)  and  base  the  solution  on  the  prescribed  displacements. 

3.  Effect  of  Damping 

In  lumped  spring-mass  models  subjected  to  transcendental  forcing  such  as  sinusoidal,  which 
occurs  within  many  electrical  con^onents  as  in  a  gas  turbine  engine,  hysteretic  damping  can  be 
included  by  considering  the  spring  stiffaesses  to  be  of  the  complex  mode  (Inman  and  Andry 
1980),  Le.,  K*  =Kr  +  iKi,  where  the  superscript  *  represents  a  complex  number,  the  subscript  r  is 
the  real  part  of  the  conq)lex  number,  and  the  subscript  /  is  the  unaginary  part.  This  representation 
is  a  mathematical  convenience  rather  than  a  physical  depiction  as  are  con^jlex  numbers 
themselves.  In  fact,  hysteretic  darqping  is  energy  dissipation,  which  is  independent  of  firequency 
and  is  proportional  to  the  square  of  the  displacement.  Hence,  hysteretic  damping  is  a  much  better 
predictor  of  the  ability  of  metals  and  polymers  to  dissipate  energy  than  is  linear  damping.  Shock 
response  spectra  are  usually  specified  with  dancing  as  a  parameter.  If  we  are  to  use  (14),  we 
should  select  the  value  of  A(o),)  for  the  appropriate  damping  value.  If  dancing  is  measured  by  the 
loss  factor  given  as  r|,  then  (14)  can  be  expressed  as 


(18) 


where  Ki  can  be  given  as. 
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(19) 


We  are  now  left  with  the  question  of  how  to  find  rii  for  each  respective  mode  of  the 
vibrational  spectrum.  In  order  to  adequately  address  this  question,  we  should  first  consider  the 
simple  case  of  a  single-degree-of-fi'eedom  system  subjected  to  sinusoidal  forcing.  We  are  then  led 
to  the  well-known  equation  of  motion  as, 

mx*  •>rk*x*  =/*,  (20) 

where 

f(t)=  Re\fe‘~].  (21) 

x(f)=  (22) 

If  we  let  the  loss  factor  to  be  defined  as 

11  =  (23) 

then  (20)  simply  becomes, 

MX*  +Kj,(^  +  ir\)x*  =f\  (24) 

Hence,  we  can  now  proceed  to  an  n-degree-of-freedom  system  such  that 

MX  *  +  iKiX  *  +  K^X  ’  =  (25) 

such  that 
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F*=Ky, 


(26) 


and 


Y{t)  =  Re\Y*^‘"  \ 


(27) 


If  we  desire  a  generalization  of  (4)  such  that  (Inman  and  Andry  1980) 


(28) 


where  is  the  real  modal  matrix  derived  from  the  undamped  eigenvalue  problem,  then  we  are  led 
to  the  following  i^*  degree-of-freedom  relationship  of  (25)  by  using  (2)  and  (3),  i.e., 


MM  +KM  + 


(29) 


Hence,  if  we  con^are  (29)  and  (24),  a  modal  loss  factor  can  be  suitably  distinguished  as 


Tl.= 


(30) 


such  that  we  are  assuming  that  the  dancing  is  small  so  that  no  significant  errors  will  arise  from 
using  the  undanq)ed  mode  shapes.  We  can  write  (30)  in  another  form,  such  as 


n,  =  ^ - ,  (31) 

is:  I 


where  we  have  the  imaginary  and  con5)lex  part  of  the  con5)lex  stiffriess.  If  we  further  let 
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(32) 


in  such  a  way  that  (31)  is  transformed  into 

Tl,-  =  - ,  (33) 

then  it  becomes  clear  that  (33)  is  the  ratio  of  potential  energies,  i.e.,  the  numerator  is  the  sum  of 
the  energies  dissipated  in  each  spring  in  mode  i,  and  the  denominator  is  the  sum  of  the  energies 
stored  in  each  spring  m  mode  i.  Hence,  this  equation  should  be  able  to  be  used  to  estimate  modal 
loss  factors,  since  it  can  be  justified  merely  on  the  basis  of  natural  intuitive  estimation.  If  we  look 
at  the  relationship  between  (33)  and  (30),  we  are  formally  led  to  the  conclusion  that  (33)  is 
essentially  the  modal  superposition  using  conqjlex  stifihiesses. 

4.  Conclusion 

This  technical  treatise  has  exemplified  a  response  spectrum  method  that  is  particularly  suited 
for  spring-mass-danping  systems  excited  by  one  or  more  base  or  ground  disturbances 
characterized  by  displacements.  It  offers  an  alternative  to  the  conventional  response  spectrum 
method,  which  is  best  suited  for  a  single  prescribed  ground  acceleration.  This  methodology  has 
also  introduced  a  method  to  incorporate  hysteretic  damping  through  con5)lex  spring  stiffiiesses. 
The  response  spectrum  method  was  formulated  in  terms  of  base  displacements  rather  than  single 
base  acceleration,  which  is  characteristic  of  base  shear  calculations  for  budding  and  structural 
vibration  analysis  under  seismic  and  ground  disturbances  or  artillery  base  loads.  When  the 
formulation  is  expressed  via  base  displacements,  a  simpler  application  of  the  response  spectrum 
method  to  varying  spring-mass  systems  with  attendant  input  excitations  can  be  realized,  as 
expressed  as  the  ratio  of  energy  considerations  inherent  in  the  natural  physics  of  the  phenomena 
itself. 
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Hence,  the  energy  ratio  equation  should  be  able  to  be  used  to  estimate  modal  loss  factors, 
since  it  can  be  justified  merely  on  the  basis  of  natural  intuitive  estimation.  Furthermore,  the  use  of 
this  derived  relation  might  be  able  to  be  applied  to  handle  single  frequency  excitations  that  occur 
in  a  typical  gas  turbine  engine  for  diagnostics  involving  component  failures  or  even  routine 
maintenance.  This  treatise  has  identified  the  original,  practical  use  of  simulating  complex 
stiffiiesses  for  hysteretic  damping,  which  is  quite  different  from  traditional  representations  of  the 
usual  spring-mass-damping  phenomenology. 
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Intentionally  left  blank. 
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3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report  will 

be  used.) _ _ _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Liformation  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to  organization, 
technical  content,  format,  etc.) _ _ _ _ _ 


Organization 

CURRENT  Name  E-mail  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7 .  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the  Old 
or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  maU.) 
(DO  NOT  STAPLE) 


